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^ Abstract. A detailed study of chemical freeze-out in nucleus-nucleus collisions 

at beam energy 158^4 GeV is presented. By analyzing hadronic multiplicities 
within the statistical hadronization approach, the chemical equilibration of 
p-p, C-C, Si-Si and Pb-Pb systems is studied as a function of the number of 
f*^ ' participating nucleons in the system. Additionally, the Two Component statistical 

' hadronization model is applied to the data and is found to be able to explain the 

observed strangeness hadronic phase space under-saturation. 
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, ■ PACS numbers: 24.10.Pa, 25.75.Dw 
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CO 1. The statistical hadronization model 

(N 

The main idea of the SHM is that hadrons are emitted from regions at statistical 
equilibrium. No hypothesis is made about how statistical equilibrium is achieved; this 
can be a direct consequence of the hadronization process pp. In a single collision event, 
, there might be several clusters with different collective momenta, different overall 

charges and volumes. However, Lorentz-invariant quantities like particle multiplicities 
are independent of clusters' momenta. 

Depending on the system size, statistical analysis can be done in different 
^ , statistical ensembles, which have been studied carefully by different groups and 

methods suitable for relativistic nucleus-nucleus collisions are well established |2|. 
In this work, small systems with net baryon number less than 10 are calculated in 
the iJS'Q-canonical ensemble taking into account exact conservation of B, S and Q 
charges. For the systems with the number of participants (Np) between 10 and 100, 
' exact conservation of strangeness only is taken into account while larger systems are 

treated grand-canonically (GC). 

In spite of the appropriate ensemble, theoretical multiplicities are calculated 
within the main version of the statistical model by fitting temperature T, scaling 
volume V, and strangeness suppression factor 75 0. Additionally, one needs to 
introduce a chemical potential for all charges that are treated in a GC manner. Usually 
the baryon chemical potential is taken as a free fit parameter, while the strangeness 
chemical potential is fixed by assuming net strangeness neutrality in the system and 
the chemical potential for electric charge by assuming that Q /B in the system equals 
Z/A of the colliding nuclei. 

The overall multiplicity to be compared with the data, is calculated as the 
sum of primary multiplicity and the contribution from the decay of heavier hadrons: 
('^j) = (nj)P'''™''''y + X]fc ^'^(^ ~^ j)("-fc), where the branching ratios are taken from 
the latest issue of the Review of Particle Physics j4j . 
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Parameters 


SHM{7s) 


SHM(7s) 


SHM(TC) 


SHM(7s) 


SHM(TC) 




p-p 158A GeV 


C-C 158A GeV 


Si-Si 15871 GeV 


T (MeV) 


177.3±5.2 


165.7±4.1 


170±10 


163.0±4.7 


162.0±7.6 


MS (McV) 




248.1±12.5 




245.5±11.0 


234.4±22.5 


7s 


0.445±0.020 


0.575±0.042 


1.0 (fixed) 


0.664±0.050 


1.0 (fixed) 


V 


0.128±0.005 


0.84±0.05 


0.23±0.03 


2.10±0.13 


0.91±0.11 








6±0.4 




11.4±1.8 


xVdof 


13.0/7 


3.4/4 


5.8/5 


7.6/4 


1.0/4 



Table 1. Summary of fitted parameters {V'=VT^ exp[-0.7GeV/T]) at top SPS 
beam energy in the framework of the SHM(75) model and SHM(TC). 

2. Experimental data set and analysis results 

The experimental data consists of measurements made by the NA49 collaboration in 
central p-p, C-C, Si-Si and Pb-Pb collisions at beam momenta 158^1 GeV [H El 13 
ISl El El d ■ The analysis has been carried out by looking for the minimum of the 

X = X^i — — ■ The fitted parameters within the main scheme SHM(75) are 
shown in Tabled A major result of these fits is that 75 is monotonically increasing 
as a function of the number of participating nucleons, and significantly smaller than 
1 in all cases (see FigQ, thus strangeness seems to be under-saturated with respect 
to a completely chemically equilibrated hadron gas. This confirms previous findings 
nails, 14, 15 . 

2.1. Superposition of NN collisions with an equilibrated fireball 

In the Two Component model (SHM(TC)), first introduced in ^21) the observed 
hadron production is taken as the superposition of two components: one originated 
from a large fireball at complete chemical equilibrium at freeze-out, with 75 = 1, and 
another component from single nucleon-nucleon collisions. Since it is known that in 
NN collisions strangeness is strongly suppressed, the idea is to ascribe the observed 
under-saturation of strangeness in heavy ion collisions to the NN component. 

With the simplifying assumption of disregarding subsequent inelastic collisions 
of particles produced in those primary NN collisions, the overall hadron multiplicity 
can be written then as (rij) = {Nc){nj) + {n'j)v, where {nj)NN is the average 
multiplicity of the j^^ hadron in a single NN collision, (Nc) is the mean number of 
single NN collisions and {nj)v is the average multiplicity of hadrons emitted from the 
equilibrated fireball. To estimate {nj)^^ in up and nn collisions, the parameters (see 
Table^ of the statistical model determined in pp are retained and the initial quantum 
numbers are changed accordingly. Theoretical multiplicities have been calculated in 
the canonical ensemble, which is described in detail in ref. |lt)| . 

This model was seen to be able to reproduce the experimental particle 
multiplicities measured in the Pb-Pb collisions at 158^1 GeV ^1 EI- For the Si- 
Si system, T, V, iib of the central fireball and (Nc) were fitted using the S-canonical 
ensemble in the central fireball. The fit quality is significantly improved compared to 
the main version of the statistical model, if the number of "single" NN collisions is 
about 11 with a 16% uncertainty. 

The central fireball produced in the C-C system needs to be analyzed in the BSQ- 
canonical ensemble taking into account the actual proton- neutron configurations. The 
Two Component model can be used to describe the C-C system as well, if the total 
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baryon number in the central fireball is B = Np + ~ 4. Even though fit quality 
is slightly worse than with the main version of the statistical model, experimental 
multiplicities are well reproduced with the SHM(TC) also. 
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Figure 1. The strangeness non-equilibrium parameter 73 and the corresponding 
\s factor as a function of the number of participating nucleons in the collision 
system (Pb-Pb from 1121 '). Both lines are of the functional form /(Np) = 
A — e~^^P and are plotted to guide the eye. Whenever the resulting ^ jdof > 1 
in the fit, errors have been re-scaled by factor \p^~[dof , for details see l4l 112). 
From left to right: p-p, C-C, Si-Si and Pb-Pb. 



3. System size dependence and conclusions 

Our fit results show non-trivial system size dependence of chemical equilibration and 
characteristic thermal parameters of the source in ultra-relativistic heavy ion collisions. 
The chemical freeze-out of all the different systems, C-C, Si-Si and Pb-Pb with beam 
momenta 158 A GeV seem to happen at similar chemical state, all of them at /is «250 
MeV. Such weak system size dependence of the baryon chemical potential has been 
already reported earlier |15| . 

Systems with few participating nucleons seem to decouple at slightly higher 
temperature than heavy systems, but in general the chemical freeze-out temperature 
as well as the baryon chemical potential are determined mostly by the beam energy, 
not by the number of participants, and thus C-C, Si-Si and Pb-Pb with the same 
beam momenta, seem not to follow the chemical freeze-out curve (see Fig. but 
show more complex system size dependence, an interplay of the initial beam energy 
and the number of participants in the system. 

The chemical equilibration of strangeness, on the other hand, seems to be strongly 
dependent on the number of participants. Going from small to large systems, the 75 
parameter increases monotonically from 0.45 in the p-p to 0.9 in the Pb-Pb with 
the same beam momenta, and thus strangeness seems to be out of equilibrium in all 
collision systems studied at SPS. The Wroblewski variable As = 2(ss)/((uii) -I- (dd)), 
the estimated ratio of newly produced strange quarks to u and d quarks at primary 
hadron level, features very similar suppression in strangeness production, see Fig. Q 
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Figure 2. Chemical freeze-out points in the [/^s — T] plane in various heavy ion 
collisions. Full round dots flgj refer to the Au-Au at 11.6 and Pb-Pb collisions 
at 40, 80, 158A GeV, whilst the square dot has been obtained by applying the 
statistical model to the preliminary 7r+, tt~ , K~^, K~ and A^p multiplicities 1181 
in the Au-Au collisions at \/S]VJV ~ ^^0 GeV. The hollow round dots refer to the 
C-C and Si-Si collisions at 158A GeV. Whenever the resulting jdoj > 1 in the 
fit, errors have been re-scaled by factor \/ I do^ ^ for details see HI 1121 . 



References 



[1] F. Becattini, hep-ph/6410403 (2004). 

[2] A. Keriinen et al, Phys. Rev. C 65 (2002) 044901 and references therein. 

[3] P. Koch, B. MuUer and J. Rafelski, Phys. Rep. 142 (1986) 167. 

[4] K. Hagiwara et al, Phys. Rev. D 66 (2002) 010001. 

[5] J. Bachler et al, NA49 Coll., Nucl. Phys. A 661 (1999) 45 

A. Billmeier, NA49 Coll., Ph. D. Thesis, Univ. of Frankfurt (Germany) 2001. 

[6] C. Alt et al., nucl-cx/0406031 and private communication with M. Gazdzicki. 

[7] S. V. Afanasicv et al, NA49 Coll., Phys. Lett. B 491 (2000) 59, J. Phys. G 27 (2001) 367, Phys. 
Lett. B 538 (2002) 275, Phys. Rev. C 66 (2002) 054902. 

[8] A. Mischke et al Nucl. Phys. A 715 (2003) 453. 

[9] T. Anticic et al, Phys. Rev. Lett. 93 (2004) 022302. 

[10] V. Friese, NA49 Coll., Nucl. Phys. A 698 (2002) 487. 

[11] C. Alt et al, NA49 Coll., nucl-ex/0409004 

[12] F. Becattini et al, Phys. Rev. C 69 (2004) 024905. 

[13] F. Becattini et al, Phys. Rev. C 64 (2001) 024901. 

[14] J. Cleymans, B. Kaempfer and S. Wheaton, Phys. Rev. C 65 (2002) 027901. 

[15] J. Cleymans et al (2002) hep-ph/0212335l 

[16] F. Becattini, G. Passaleva, Eur. Phys. J. C 23 (2002) 551. 

[17] J. Manninen et al, nucl-th/0405015 

[18] D. Rohrich, talk given in Strangeness in quark matter 2004 meeting. 



